Criteria for shear banding in time-dependent flows of complex fluids 
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Within a highly generalised theoretical framework for the flow properties of complex fluids, we 
study the onset of shear banding in the three most common time-dependent experimental protocols: 
step stress, step strain and shear startup. By means of a linear stability analysis we derive a fluid- 
universal criterion for the onset of banding, separately for each protocol, that depends only on the 
shape of the experimentally measured time-dependent rheological response function, independent of 
the constitutive law and internal state variables of the particular fluid in question. Our predictions 
thus have the same status, in these time-dependent flows, as the widely known criterion for banding 
in steady state (of negatively sloping shear stress vs. shear rate). We support them with simulations 
of the rolie-poly model of polymeric fluids, the soft glassy rheology model, and a fluidity model. 

PACS numbers: 61.25.he,83.50.Ax,62.20.F,83.10.y,83.60.Wc 



Many complex fluids show shear banding [1,] , in which 
an initially homogeneous sample of fluid separates into 
layers of differing viscosity under an applied shear flow. 
Examples include wormlike micellar surfactants 
polymers soft glassy materials [3, H|, and (possi- 
bly) biologically active fluids @ . At a fundamental level 
shear banding can be viewed as a non-equilibrium (flow- 
induced) phase transition, or equivalently as a hydrody- 
namic instability of viscoelastic origin. In practical terms 
it drastically alters the rheology (flow response) of these 
materials and thus impacts on widespread industrial ap- 
plications: in plastics, foodstuffs, well-bore fluids, etc. 

In steady state, the criterion for shear banding is (usu- 
ally 0) that the underlying homogeneous constitutive 
curve of shear stress vs. shear rate has negative slope, 
dTi/dj < 0. Beyond studies in steady state, accu- 
mulating data in polymers @-15|, surfactan ts jliWlij l. 
soft glassy materials flil - Eil ]. and simulations [23M30j ] re- 
veals that shear bands often also arise during the time- 
dependent process whereby a steady shear flow is estab- 
lished from an initial rest-state. These "transient" bands 
can be sufficiently long lived to represent the ultimate 
flow response of the material for practical purposes even 
if the constitutive curve is monotonic dY./dj > 0. 

In view of these widespread observations, crucially 
lacking is any known criterion for the onset of banding in 
these time-dependent flow protocols, of the same status 
as given above in steady state: independent of the consti- 
tutive model and internal state variables of the particu- 
lar fluid in question, and depending only on the shape of 
the experimentally measured rheological response func- 
tion. Here we provide such criteria, separately for each 
of the three most common protocols: step stress, step 
strain and shear startup. We do so within a highly gener- 
alised framework that includes many widely used models 
of polymeric fluids (polymers and wormlike surfactants) 
and soft glassy materials (foams, dense emulsions and 
colloids). Our aim is thereby to develop a unified un- 
derstanding of existing and future observations of time- 



dependent banding; and to facilitate the design of proto- 
cols that optimally enhance or mitigate it. 

Consider a sample of fluid sandwiched between par- 
allel plates at y — {0,L}, sheared for times t > in 
one of the protocols defined below. We decompose the 
shear stress in any fluid element into a purely Newto- 
nian part of viscosity rj, and a viscoelastic part due to 
the polymeric or soft glassy degrees of freedom: S(t) = 
Ga(y,t) + ryy(y,t). Transnational invariance is assumed 
in the flow direction x and vorticity direction z. Het- 
erogeneity is allowed in the flow gradient direction y, 
though force balance (neglecting inertia) demands uni- 
form E. The viscoelastic shear stress Ga(y, t), written as 
a modulus G times a dimensionless conformation variable 
cr, obeys dynamics in which loading by shear competes 
with intrinsic stress relaxation. In the simplest descrip- 
tion dt<J = 7 — -cr, though in reality this shear compo- 
nent cr = a X y also couples to normal stress components 
Vxx , &yy ■ ■ ■ |3l| • In models of soft glasses it also couples 
to (at least one) variable a characterising the sample's 
fluidity [27, 32, 33 1. Collecting all such variables into a 
vector s = (cr, a xx ■ ■ ■ ) T and defining the projection vec- 
tor p = (1,0, • • • ), we write in the most general case: 



E(i) = Gp.s(y,t) 
d t s(y,t) = Q(s,7). 



(i) 

(2) 



For specific choices of s and Q, Eqns. (jTl [2j) encompass 
the predictions in shear of most commonly used rheo- 
logical models. Within this general framework we will 
derive fluid-universal criteria for the onset of banding, 
separately for step stress, step strain and shear startup 
protocols. We support these with simulations of the 
rolie-poly model |34ll of polymeric fluids; a fluidity model 
of soft glasses [23]; and the soft glassy rheology (SCR) 
model [27], [35[. (All these are defined in the Suppl. Mat.) 

Terms of the form DV 2 s must also be added to Eqn.[2] 
to prohibit flow heterogeneity at scales smaller than the 
fluid microstructure I — VDt, and to ensure unique se- 
lection of any shear banded state [36| . In any linearised 



FIG. 1: a) Creep curve in the non-stretch rolie-poly model for 
ft — 0.8, n = 10~ 4 , Eo = 0.7. Dashed: linearly unstable; solid: 
linearly stable, b) Corresponding degree of shear banding, c) 
Flow profiles at times marked by circles in a) for £„ = 0.1(5n,o, 
I — 10 -2 with increasing v(y — 1) for increasing time t. 



equation, e.g. (JIJ) below, these can be neglected for the 
fluctuations of interest here with wavelength ^S> I. 

We use units in which the modulus G = 1; the rheome- 
ter gap L — 1; and (as defined in the Suppl. Mat.) the 
underlying fluid relaxation timescale td = 1 (rolie-poly), 
or to = 1 (SGR and fluidity models). 

Step stress protocol — Consider first an experiment in 
which a step stress of amplitude Eo is applied to the 
sample at time t — 0. We express the response to this 
load as a spatially uniform "base state" of homogeneous 
creep, plus any (initially) small heterogeneous departure 
from it: j(y,t) = j (t) + J2 n S in cos(rar?//L) exp(wt), 
s(y,t) = s (t) + J2n Ss nCos{mry/L)exp(ujt). (Note, 
then, that we are studying the time-evolution of pertur- 
bations to a base state that is itself time-dependent 37].) 
To investigate whether at any time t the system is linearly 
unstable to the onset of banding, we expand Eqns. (JTJJSJ) 
in powers of the heterogeneous perturbation's amplitude 
and neglect terms beyond first order. 

The dynamics of the homogeneous base state is speci- 
fied by the zeroth order terms as follows: 



d t (d t s ) = [M — Gqp/rj] ■ d t s Q 



(3) 



with M = <9 s Q|s ,7o and q — d^Q\ 8o ^ . This equa- 
tion must be solved subject to the initial condition 

9tSo|t=o = Q(so|t=0)7o(* = 0) = £ /»7). Tne dynam- 
ics of the heterogeneous fluctuations is specified by the 
first order terms according to 



d t Ss n = [M — Gqp/rj] ■ Ss r , 



(4) 



This must be solved subject to source terms specifying 
the seeding of any heterogeneity, whether due to (i) sam- 
ple preparation protocol, (ii) slight flow device curvature, 
(iii) mechanical or (iv) thermal noise. Here we consider 
(i), modelled as an initial condition Ss n (0) = e n N n with 
each e n a small amplitude and the entries of N n chosen 
randomly from a flat distribution of mean and width I . 

Eqns. (J3j HJ) reveal that the heterogeneous fluctuations 
Ss n obey, in this step stress protocol, the same dynam- 
ical equation as the time derivative dtSQ of the homoge- 
neous base state. Due to their different initial conditions, 




FIG. 2: a) Creep curves of the SGR model for stress values 
So/S y = 1.005, 1.010- ■• 1.080 (curves upwards), b) Corre- 
sponding degree of shear banding, c) Normalised velocity 
profiles for the circles in a), x = 0.3, w = 0.05, n = 50, m = 
10000. Initial sample age t„ = 10 3 [1 + ecos(27rj/)], e = 0.1. 



however, dtSo and 6 Sfi circ not guaranteed to evolve col- 
inearly. Numerically, though, we find they always do be- 
come colinear after a very short transient. Accordingly 
any component of Ss n grows whenever its counterpart 
in So obeys dfso/dtSo > 0. Combined with the linear- 
ity of (JTJ this shows that banding must develop when- 
ever the underlying homogeneous creep function 70 (i) 
obeys d^'yo/dt'jo > 0. This important result, valid for 
any model of the highly general form (JTJ |5J), is consistent 
with numerous recent experiments in polymeric fluids 0- 
13, EMI El and soft glassy materials {2lll2^. 

We support our prediction with numerical results for 
the rolie-poly model 34( of polymeric fluids in Fig. JTJ 
with parameters for which the constitutive curve £(7) 
is monotonic and the steady state unhanded. A homo- 
geneous creep function 70 (t) is shown in Fig. [1^, with 
the regime of linear instability to the onset of band- 
ing, <9j7o/<9 t 7o > 0, dashed. The corresponding degree 
of banding as measured in a full nonlinear simulation 
(FigJTJj) is consistent with this prediction; and snapshot 
velocity profiles confirm obvious banding as the system 
starts flowing under the applied load (Fig. JTJ:). 

We also simulated the spatially aware soft glassy rhe- 
ology model 27, [35| in its glass phase, for which the 
steady state constitutive curve shows a yield stress £ y 



with monotonic increase beyond: £ = £ y + aj m with 
to > 0. For an applied stress just above £ y we see an 
early time regime of slow creep in which the shear rate 
progressively decreases (Fig. J5^). This then gives way 
to a steady flowing state via a regime of upward slope 
dtj > and curvature df"j > in which banding arises 
(Fig. J2b,c), consistent with our prediction above. Ho- 
mogeneous flow is recovered in steady state, though for 
stress values only marginally above yield the recovery 
time is likely to be unattainable experimentally: band- 
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FIG. 3: Left: stress vs. strain for a fast ramp in the rolie- 
poly model, f} = 0.0, tr = 10 _4 ,7? = 10~ 5 . Right: relaxation 
post-ramp; unstable region dashed. Upper curve: appreciable 
chain stretch, trj oo. Lower: negligible stretch, TR-y = 0.1. 



ing will then represent the material's ultimate observable 
flow response in practice. 

Step strain protocol — Next we consider a rapid shear 
strain ramp 70 = jot applied for times < t < t* , af- 
ter which the strain is held constant at 7q = jot* ■ The 
limit 70 — > 00 at fixed 7g gives a true step. As above we 
analyze the response to this deformation by first (artifi- 
cially) restricting the flow to remain homogeneous, then 
performing a linear stability analysis for the dynamics of 
small shear banding fluctuations about this base state. 

Assuming decomposition of the viscoelastic dynamics 
@ into additive loading and relaxation parts, jS and 
R/t, the base state obeys during the ramp 



d l0 s = S(s ) - ^R(a ) w S{s Q 



(5) 



dominated by loading for j Q r 3> 1. Post-ramp it relaxes 
back to equilibrium according to dtSo = — -R(sq). 

To investigate whether this relaxation can indeed oc- 
cur in a purely homogeneous way, or whether the sam- 
ple instead transiently shear bands during it (assum- 
ing t* so short that no banding arose during the ramp 
itself), we write j(y,t) — J2 n ^7" cos(n7ry/L) exp(wt), 
s(y, t) = so(t) + J2 n cos(niry / L) exp(wi) and expand 
to linear order to find that the heterogeneous fluctuations 
evolve post-ramp according to 



d t Ss r , 



£S(s )p-±d s R\. 



■5s n w -£S(s )p-6s n 
(6) 

for small solvent viscosity 77 -C Gt. 

Denoting by Sq the system's state instantaneously as 
the ramp ends, we see from (|5|) that Sq relates to the evo- 
lution of the base state just before the ramp ended via 
0~so|t=t*- = S( s o)- Because Sq is continuous at t* , Sq 
also relates (by [6j) to the dynamics of the fluctuations im- 
mediately post-ramp via d t 6s n \ t=t ,+ — — j-S(so)p.8s n . 
Combining these, taking the first component, and ap- 
pealing to the linearity of force balance, we find 



dt5jn\t=r+ = --^^-\ t =t*-5in- 

7] (77 



(7) 



This result tells us that, as the system commences its 
relaxation back to equilibrium, it is unstable to banding 



if the stress had been decreasing with strain just prior to 
the ramp ending. The prediction applies for any rheologi- 
cal model with additive loading and relaxation dynamics, 
and is consistent with early intuition (39j . 

Numerical results for the rolie-poly model (Fig. J3J) sup- 
port this. As discussed in the Suppl. Mat., this model 
has two intrinsic relaxation times: that on which a poly- 
mer chain escapes its tube of entanglements, td, and that 
on which chain stretch relaxes, tr -C td . The lower curve 
in §3§ is for a ramp rate 7TD 3> 1, 7Tr <C 1. This corre- 
sponds to the 7 — > 00 limit of Eqn. [TS] with its nonlinear 
loading terms: during the ramp the material behaves as 
a nonlinear elastic solid with a maximum of stress vs. 
strain. If the total applied strain 7* exceeds the location 
of this maximum, the system is left unstable to band- 
ing immediately post-ramp. The upper curve in ((3]) is 
for a much faster ramp 7TR 3> 1. This corresponds to 
the 7 — > 00 limit of Eqns. 1161 which have linear loading 
terms. Here the material behaves as a linear elastic solid 
during the ramp and is stable immediately after it. 

The upper curve reveals further polymer physics. Re- 
laxation of chain stretch on the short timescale tr, post- 
ramp restores a state as if no stretch had arisen in the first 
place: the upper curve rejoins the lower, both are unsta- 
ble to banding and only finally decay on the timescale td- 
These results are consistent with experiments [l3- T^, 40l | 
and numerics 24,[26| showing that banding can either oc- 
cur straight after a step strain, or following an induction 
period. In extensional equivalent, it might also underlie 
the physics of delayed necking instability [ill Ell ■ 

The fluidity and soft glassy rheology models have lin- 
early increasing stress during a fast ramp so are predicted 
by Eqn. ([7|) to be stable against banding after it. 

Shear startup - We turn finally to a shear startup ex- 
periment in which a constant strain rate 70 is applied 
for all times t > 0, giving strain 70 = jot. As always 
we study a base state given by the homogeneous solu- 
tion {jo,so(t)} of Eqns. ((TJ ^ in this protocol, and the 
dynamics of small banding fluctuations about it. 

Consider first a startup run in which the sample is 
artificially constrained to remain homogeneous until it 
attains steady state t — > 00. In this limit (with the con- 
straint now removed) the criterion for instability is known 
to be simply dj T.Q < 0. To generalise this to any finite 
time during startup we start by taking 9^, | 7o [(*)/7o]j 
where (*) denotes Eqn. [21 to find that the base state 
obeys 

9jo a oUo = P ■ M ~ X ■ [^to^oIto + TcAtAoSo - q] , (8) 

with M, q defined as above. Heterogeneous fluctuations 
about this base state obey Eqn. [4] the largest eigenvalue 
of which crosses zero when | — ^qp + M\ = 0. Com- 
bined with ((5J) this gives, after some manipulation, the 
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FIG. 4: Rolie-poly model with a non-monotonic constitutive 
curve (/3 = 0.4, tr = 0.0,77 = 10 -4 ). Any horizontal slice 
across this plane corresponds to a single shear startup run. a) 
Large circles: full onset criterion Right dashed line: vis- 
cous criterion (|10[l . Left dashed: elastic criterion <|12[). Dot- 
ted: stress overshoot d 70 £o = 0. Small circles: elastic plus 
viscous terms (|12|) + (|10|) . b) Circles as in a), and contour lines 
of equal S'jn = 7ol0 M for integer M, from integrating Eqn. [4] 
with initial perturbation e n — 10~ 2 . First contour: M — —2. 



condition for banding instability during startup as 

|M|{9 70 £ | 70 -P-M 1 ■ [a 7o s| 70 +7o^ 9 7o s]} < 0. 

(9) 

Note this contains both "viscous" derivatives with re- 
spect to strain rate <9 7o (as in the steady state criterion) 
and now also "elastic" ones <9 7o with respect to strain. 

To explore the physical meaning of this result we con- 
sider now a whole series of startup runs: each at a differ- 
ent shear rate 70 , and with the evolution of the base state 
So monitored as a function of strain 70 (or equivalently 
time t) in each. Two limiting cases are instructive. 

First is the steady state limit t — > 00 of any individual 
run. Here the accumulated strain 70 becomes irrelevant 
and (O reduces to the condition for steady state bands: 



|M| 9 70 S | 70 < 0, 



(10) 



with just the viscous derivative <9 7o . (Note that \M\ > 0.) 

The second limit concerns startup runs performed at 
very high rate 70 — > 00. If the stress startup curve attains 
a unique limiting function £(70) independent of 70 in this 
regime (and likewise for all other components of s), the 
material effectively then acts as a nonlinear elastic solid. 
Only elastic derivatives <9 7o then act in © to give: 

-\M\p-M- 1 -d 10 s\ % < 0, (11) 
-trMS 7o S | 7 o+7o^ Eo| 7o < 0, (2D) (12) 

in which trM < 0. Here (fl~2"|) specialises ([Til) to the case 
of just two degrees of freedom in s as relevant, e.g., to 
the non-stretch rolie-poly model for which s = (a, <J yy ) T . 

These predictions are illustrated in Fig. U by numerical 
results for the non-stretch rolie-poly model with param- 
eters for which the constitutive curve is non-monotonic 
and the steady state banded. Any horizontal slice across 




FIG. 5: Shear startup in a) rolie-poly model with a mono- 
tonic constitutive curve. /3 = 1.0, tr = 0.0,7? = 10 -4 . Large 
circles: full onset criterion ([9|. Small circles: elastic plus vis- 
cous terms (|12I) + H10I) . Solid lines: contours of 5j n = 7ol0 M , 
from integrating Eqn. [4] with initial perturbation e n = 10~ 2 . 
First contour: M — —2. b) Fluidity model for 77 = 0.05, ini- 
tial age t w = 10 4 . Circles: onset criterion ((9|. Dashed line: 
<9 70 £o = 0. e n as for a) (but with perturbation added to a 
only). In a) and b) squares delimit region of shear banding, 
(7/1 — > 0.05, in nonlinear simulation with I — 10 -2 . 



this (jo,t) plane corresponds to a single startup run. 
"Elastic" instability (fl~2"|) is apparent along any such slice 
at high strain rate. The first term of (|12l) . taken alone, 
predicts onset of banding at the overshoot in the startup 
curve So (70); the second adds a correction to cause onset 
just before overshoot, as seen numerically. Steady state 
"viscous" instability (fTTJ|) is seen at the far right t —> 00. 
Between these limits the two instabilities combine and 
in principle one needs the full criterion Q , of which the 
cross terms <9 7o 9 7o would be hard to measure experimen- 



tally. Pleasingly, however, simply combining criteria (|10p 
and ([T2|) (i.e., setting the sum of their LHS's negative) 
without the cross terms of (JH) captures onset in most 
regimes, except for a small region at the lowest 70. For 
a monotonic flow curve (Fig. [5^,), steady state instabil- 
ity is absent but a diagonal patch of transient elastic-like 
instability remains. 

These results are consistent with experiments [1, [l^, 
lil . HI and numerics 2^, 24 , 26 , 13 for polymeric systems 
showing transient bands after stress overshoot, which 
may or may not persist to steady state depending on 
the slope of the ultimate flow curve. 

The SGR and fluidity models do not attain a limiting 
startup curve £0(70) but retain dependence on 70 even as 
70 — > 00. Separation into pure viscous and elastic insta- 
bilities is then not possible, though the steady state con- 
dition for banding <9 7o £o < remains. Analytics for the 
fluidity model (which has <9 70 So > and no steady state 
banding) nonetheless predict a startup instability onset 
immediately beyond stress overshoot <9 7o £o < 0, consis- 
tent with experiment 3, 2(| . Numerics further show the 
bands to be very long lived at low j , Fig. [Sp 44]. The 
SGR H3 and STZ [IE [HI models (not shown) likewise 
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show long-lived bands onset at overshoot. It remains 
unclear whether this onset criterion 9 7o So < holds in 
any model (such as these) for which the stress depends 
strongly on strain but only logarithmically on strain rate. 

Conclusion - We have given fluid-universal criteria for 
shear banding in time-dependent flows of complex flu- 
ids. In a step stress protocol, banding is predicted in 
any regime where the shear rate of the underlying creep 
response obeys dfjo/dtjo > 0. For a fast strain ramp 
("step strain") protocol, banding instability is predicted 
immediately post-ramp if the stress had been decreas- 
ing with strain just prior to ending the ramp. In shear 
startup we predict separate "elastic" and "viscous" in- 
stabilities, and a "mixed" region between these, for a 
broad class of fluids that attain a limiting nonlinear stress 
startup curve £(70) in fast flows 70 — > 00. In contrast the 
startup curve of soft glassy models retains dependence on 
70 even as 70 — > 00 and these fluids show a startup in- 
stability that is neither purely elastic nor viscous, onset 
just after stress overshoot. We hope these predictions will 
help unify our understanding of widespread data concern- 
ing shear banding in time-dependent flows, and stimulate 
further experiments to test the picture suggested here. 

The authors thank Stephen Agimelen, Mike Cates, 
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Smeeton and Peter Sollich for interesting discussions; and 
the UK's EPSRC (EP/E5336X/1) for funding. 
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I 

SUPPLEMENTARY MATERIAL 



Here we detail the governing equations of the models for which numerical results are presented in the main text. 
In each case we adopt the continuity equation for incompressible flow: 

V-tt = 0, (13) 

in which u is the fluid velocity field. We also impose force balance at zero Reynolds number: 

V • (Go- + 2rjD) = 0, (14) 

in which D = \{K + K T ) with K = Vix. Eqn. [14] corresponds to Eqn. [2] in the main text, with incompressibility 
(Eqn. 1131) automatically ensured for any flow of the form considered throughout this work, with spatial variations 
only in the flow gradient direction. The viscoelastic dynamics of cr are specified for each model separately as follows. 



Rolie-poly model of polymeric fluids 



The rolie-poly model [34j of polymers and wormlike micellar surfactants has 



da 

~dt 



K a ct K 1 — (o- - I) - — (1 - A) \tr + (3A- 25 (<t - 1)1 . 

TD T R 



(15) 



in which A = y3/T for molecular trace T = <r xx + <r yy + a zz . In this equation td is the characteristic time for a 
polymer chain to reptate out of its tube of entanglements and tr is the timescale on which chain stretch relaxes. 
Following [2H we set S = —1/2 throughout. In planar shear flow for which u — jyx, Eqn. fp~5]) has componentwise 
dynamics 



d t cr 

dtVyy 



-—a -—(l-A)(l + 0A)a, 

TD Tr 

1 , ,2 



mi 



1 



TR 



(1 - A) [<Jyy + PA{ayy - 1)] , 



TD 

270- (T - 3) - —(1 - A) [T + (3A(T - 3)] 

TD TR 



These equations fit the generalised framework in the main text with a three dimensional s — (a 7 a y y,T) T . 
In the limit of fast chain stretch relaxation tr — s- Eqn. (jT5|) reduces to 



dcr 

~dt ~ 

with componentwise dynamics 



K ■ a + cr K 

d t a = 7 
2 



1 



TD 



-tT{K.cr) [tr + p(tr-I)] 



(16) 



(17) 



1 

— ct, 
td 



TD 



(18) 



7 



and T = 3. These fit the framework in the main text with a two dimensional s = (<r, <J yy ) T . 

As noted in the main text, a diffusive term of the form £>V 2 <x must be added to Eqn. [15] in any spatially resolved 
simulation of the shear banded profiles. 



A fluidity model of soft glasses 



A simple fluidity model of soft glassy rheology has 



^ = K * + <tK t --(a -I), (19) 
at t 



with the stress relaxation time r ascribed its own dynamical evolution according to 



? = 1 (20) 

dt Tq + 1/V2tr£>.£> 

In a planar shear flow u — jyx this model has componentwise dynamics 

d t a = 7 - -a, 

T 

dtT = 1 1 —FrT7- (21) 

ro + l/|7l 1 ' 

Transforming r — > a = 1/r this fits the framework in the main text for s = (c, a) with a fluidity variable a. Spatially 
diffusive terms must again be included in any simulation of the shear banded profiles. In this case model we add these 
only to the fluidity variable r following 27] . 



The soft glassy rheology model 

The soft glassy rheology model in its original (homogeneous) form [35| considers an ensemble of elements undergoing 
activated hopping over local energy barriers E, governed by a noise temperature x. Each element is assigned a local 
strain I which, between hops, obeys I = 7 where 7 is the externally applied shear rate. The characteristic hop time of 
any element is 

[E-±kl 2 \ , 
T = r exp 2 . (22) 



After any hop, an element resets its local strain I — > and chooses a new yield energy E at random from a distribution 
p(E) — exp(-E). The evolution of the ensemble's distribution P(E,l,t) of yield energies and local strains 13511 can 
be cast exactly in terms of the corresponding dynamics of infinitely many stress- and fluidity-like moments [45j . 



dl I dE—P(E,l,t), (23) 







T 



to give 



^yn 1 

9tP m ,n = 7mim-l.fl H Prn+l.n ~ Prn,n+1 + Z — \ T-Po,lS m ,0 (24) 

X ' 1 + n/x 



for m, n = 0, 1, 2 • • • . Listing P m , n in the state vector s, and noting that P\$ — o (i.e. the macroscopic stress is the 
average of the local ones) we recognise this to fit the framework in the main text with infinite dimensional s. 

As described in Ref. [27l |43| , this model is rendered spatially aware by means of a numerical simulation in which 
we evolve j = 1 ■ • • m SGR elements on each of i = 1 • • • n streamlines, corresponding to y = • • • 1 with periodic 
boundary conditions. The stress on streamline i is <Ji — (1/m) Uj. At any timestep a waiting time Monte Carlo 
algorithm chooses stochastically the next element to jump. Supposing the jump occurs at element ij when its local 
strain is I = I, force balance is then imposed by updating all elements on the same streamline as I —> l+ljm. Diffusive 
coupling of the dynamics on neighbouring streamlines is included by further adjusting the strain of three randomly 
chosen elements on each adjacent streamline i ± 1 by £w(—l , +2, — 1). 



